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I 

The ! resonance^ below the inelastic (n=£) threshold of elec- 


tron-hydrogen scattering are examined by the projection operator 
technique of Feshbach. This technique converts the problem to an 
eigenvalue problem for a projected Hamiltonian in which the auto- 


ionization levels of H~ become true bound states. The number of levels 
of this problem is found to be infinite and cluster at the n=2 
level of H. The relative spacings of the levels are identical to 
the resonances in e-H scattering as derived by Gailitis and Damburg. 
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i nmoDUCTioN 

Recently, O'Malley and Geltman"* - have utilized the projec- 

2 

tion operator technique of Feshbach , as applied to the atomic case 

3 

by Hahn, O'Malley and Spruch , to calculate the compound auto- 
ionization states of H below the n=2 excitation threshold for e-H 
scattering. The point of the projection operator, Q, is that it 
projects out (removes ) the ground state of hydrogen from the total 
wave function, and therefore gives rise only to asymptotically vanish- 
ing terms even for a wave function that would otherwise be capable of 
describing e-H elastic scattering. Thus one can formulate an 
ordinary eigenvalue problem for the Schrildinger equation with a 
projected wave function Q : 

(i 


4 

and, following Fano , O'Malley and Geltman have shown that at 
energies t-. close to the eigenvalues °f the above problem. 


t 





( 1 . 2 ) 


the phase shift associated with the scattering problem, described 
by ~\jf , undergoes an increase of approximately 7T radians cor- 
responding to resonances in the elastic scattering. 

Although an explicit formula for /L^was given^", it involves 
coupling to the continuum part of . Thus the O'Malley-Geltman 
calculations were confined to u below the inelastic (2s) 
threshold for H (e-H scattering) and He (e-He + scattering). 

The neglect of would indeed be justifiable as agreement 


c ^ y g 1 

with other scattering calculations'^ 3 1 and experiments 0 showed X . 
However, there is a special circumstance in the case of elastic 
scattering just below the 2s threshold which makes for an infinite 
number of resonances which cluster on the 2s threshold. This 

Q 

remarkable result was deduced by Gailitis and Damburg on the 

basis of close -coupling scattering calculations above threshold 

■which were then extrapolated below threshold. The physical basis 

for this phenomenon derives from the very long range nature of 

the potential ( r~ x ) which in turn stems from the degeneracy of the 

2s and 2p levels of hydrogen (and it is only in this approximation 

that the infinity of resonances applies). On the other hand the 

actual calculations of O'Malley and Geltman^ yielded only two 

result 

resonances in each channel, which, A although it was in accord with 
previous predictions of two resonances^, leaves open the question of 
whether the higher resonances are present below threshold in the 
actual spectrum of Q 1 jf or whether they appear above the 2s threshold 
and are only brought below threshold by the un calculated 

The formal definition of the projection operator Q as it 
applies to this problem has been given by Hahn, O'Malley, and 
Sprueh 2 : 


.th 
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is the ground state of hydrogen Occupied hy the 
particle (i = 1,2). 

In the next section we shall show that the spectrum of discrftte 
levels of is indeed infinite below threshold, and we shall 

derive a formula for the relative spacings of which is 

identical to that of Gailitis and Damburg^ for 

II ANALYSIS 

We start with a close -coupling wave function for the A S state 
of e-H scattering including ls-2s-2p states of hydrogen: 



cn6, x ] 
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where is r times the given radial wave f -unction. A 

straightforward application of Q on now yields 
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(2.2a) 


where 



'or 



/P (VJ ~y) dr 

/4- 


(2.3) 


Here we note that one can add any amount of R 2S to u and simul- 
taneously subtract the same amount to R ls from v without changing \j/ . 
Thus by selecting this amount as C^- , i.e., by redefining u and v 
as given in (2.4): 


U.(C t ) 


urn) +- C- 


(2.4a) 
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n/rrj 



R Cn-> 

/A^ 
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(2.4b) 


we find, without having changed , that 



nj-(r,) R 2 l r ^ R 2 ^ r ^ Ccr3 -^ z 



(2.2b) 


Therefore without any loss of generality we can take (2.2b) as 
the general form of Q) ~yJ , and we see that it has the desired 
vanishing asymptotic form below the n=2 threshold: 

iU* QTf= 0 

f . — >> oo (2.5) 

1 L 


is in fact the close coupling wave function with the Is 
term absent. The equations which v and w satisfy may be simply 
derived from the integro-differential equations for the acattering 
in the close coupling approximation 10 by dropping all terms depending 


on u and Rxs* At large distances they reduce, for S states, to 
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These equations may he diagonalized so that two linear combinations 
of v and w satisfy two uncoupled equations. The one which contains 
the attractive r 2 potential corresponds to 

f(r) = *T(r) b ( ^ Cr) • 


The diagonalized equation is 


f b ^ (*&-') +£) FO)* o. 


(2.6a) 


(2.6b) 


( 2 . 7 ) 


( 2 . 8 ) 


In this equation 


£ 


t 


is the energy in rydhergs relative to 


the 2s threshold. Having gotten the equation in this form, we 

12 

can appeal to an argument of landau and Lifshitz to answer 
affirmatively the question of whether there are an infinite number 

"|p 

of bound states of F(r). For they have shown that any problem 
governed by the potential in (2.8) for large r and well-behaved 


for small r has a solution with an infinite number of nodes as 
— -> 0 . Thus for £1 negative and small, one can find bound 
state solutions with any arbitrary number of nodes, i.e., there 
are an infinite number of bound states with 0 

9 

This much was also clear to Damburg and Gailitis as a 
perusal of their paper shows. The new point thus far has been 
that the elimination of the Is state from the close -coupling wave 
function can be put on a more rigorous basis by a consistent 
application of the Q operator. Below we shall actually evaluate 
the eigenvalue spectrum of F. 

Letting 


6 - 



( 2 . 9 ) 


! 


one can convince oneself that the solution of (2.8) which vanishes 
at infinity is 



€ 



( 2 . 10 ) 


where 


•I(np-Gfr) 


( 2 . 11 ) 


and 



is the Hanhel function 





J and N y are the Bessel and Neumann functions respectively. 
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For large values of r 


Ffr) OC e (f 

Y—> 




are 


( 2 . 12 ) 
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For values of r < a the equations for F and a second linear 
combination of v and w, call it G = v + [(1 + \Z57)/6]w, will be 
coupled. They cannot be decoupled by any linear transformation. 

14 

However, it is certainly possible to define a non-local potential 

which will describe an uncoupled F(r) down to the origin. This 

non-local potential can in turn be approximated a local potential 

V (r) at least for values of r slightly below the value of r = a 
e 

where Eq. (2.6) is valid. Now it can safely be assumed that V e (r) 
will continue to decrease (increase in magnitude) until such a point 
r = b where the centrifugal barrier forces the potential, up to 
infinity. I.e., 


|V e (r)| ;> |V(a)| = 


b < r < a. 


(2.13) 


Then for any energy satisfying 


c « 


8.2 


(2.1k) 


the wave function F(r) will be virtually independent of e for r £ a. 
Let the logarithmic derivative of F(r) at r = a, as obtained by integ 
ration from the origin in the presence of the actual potential, be 
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defined by 


t 




(2.15) 


To repeat , for energies satisfying (2.14), t is virtually independent 

of e • 

This logarithmic derivative must be matched to that of F(r) 
for r > a in order to obtain the allowed bound states. In this 
region F(r) is given by (2.10). It is shown below that as long 
as we only consider energies consistent with (2.14), the small 
argument expansion of the Hankel function in (2.10) may be used 
in the vicinity of r = a. 

13 

To see this we consider the small argument expansion 


H 


(1) 


(z) = 


X 

sin V- 


[ e -iVTT (fcs> 

771 TTl^ 


(%z) 

Ki^y 


J • C 1 * 0 3 * (2 - 16) 


Putting in the imaginary values V = ia and z = ig 2 r yields the 



explicit expression 



2ie^ aTT sin 


1 



[to sinh (tto )]2 




[ 1+ <vki ) ] 




where 


cp = arg [r(l + ia)] . 


It is easily verified that for a given by (2.1l) and for r a, 
(2.14) guarantees that 



1 . 


Therefore, with no further restriction we have, to within a 



9 


- 13 - 


proportionality constant 


F(r) = Jr sin [a In r) 


tp], r p» a 


(2.18) 


and 


fTsj 


dF(r) 

dr 


= | cot [a In (& a) - cp] + 
r=a ^ 


( 2 . 19 ) 


The allowed hound state energies may he found hy equating this 
to t, the value from the interior, and solving for assuming, of 
course, that t is independent of e . 


i o 

e 2 = — exp 


{ : O'Xn 


ta-1 ^ 
2a J 
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+ nrr J j- 


( 2 . 20 ) 
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where 


_i 

0 £ cot X ^ TT 


Let the energy of the lowest hound state he designated hy ei . 
If satisfies (2.14), the successively higher states are given hy 
the recursion relation: 


e n+l 


e -2Tr/a 


( 2 . 21 ) 


as is easily seen from (2.20). Equation (2.2l) is identical to 

the formula for the ratio of resonant energies obtained hy Gailitis 

g 

and Damburg . This relation becomes more accurate with increasing 
n, because of the decreasing dependence of t on c , and because of the 
increasing validity of the first order expansion ( 2 . 17 )* We now look 
at the case n = 1. 

O'Malley and Geltman 1 have given for the energy of the lowest 


lying 's bound state: 



* 


» 
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ei = 0.645 ev = 0 . 0474 Rydbergs. 


( 2 . 22 ) 


If this energy is to satisfy (2.l4) ; the matching radius which we 
choose must satisfy 


a 2 

e i 


107 , 


where all lengths are given in Bohr radii. For a = 5, (2.14) is 
satisfied fairly well. However, the extent to which the asymptotic 
equations are valid at such a small radius is questionable. The 
potentials coupling in other channels could give 10$ effects at 
r = 5 quite easily. But allowing a 10$ error in the relation 
(2.2l) for n = 1 is not serious, since with ei fixed by (P. PP ) 
one obtains 


e 2 = (0.037 + 0.004) e.v. 


(2.23) 
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This res ul t corresponds to a fairly narrowly defined second 
eigenvalue; namely (relative to the ground state of hydrogen): 


£ 2 = (IO.I67 + 0.00^) e.v. 


Furthermore, since the asymptotic form of the 2s-2p equations, 
(2.6), contains no exchange terms, Eq. (2.19) must hold for triplet 
as well as singlet S states (albeit with a different ex). For 
other values of the total angular momentum one has equations similar 
to (2.6) with altered numerical coefficients of the r 2 terms 
Only in P and D states will the diagonalization yield a net at- 
tractive r 2 potential with associated coefficients a = 1.86 and 0.75 
respectively . With these changes in a, Eq. (2.19) also holds for 
states of both spins. 

In Table I, we summarize the numerical results for the first 
three resonances for S and P states. The energies are given relative 
to the ground state of hydrogen (which is appropriate for comparison 
with resonant energies in the elastic scattering of electrons from 
hydrogen). The S and P states of both spins are included and the 
lowest state in each case has been chosen to agree with the calculated 
results of O'Malley and Gelt man"*". The second states, which in the *S 
case is in accord with the original prediction of two resonances^* 7, 


f 


* 


- 17 - 

then provide a point of comparison with the 0 'Malley-Geltman result. 

It can he seen that our second states particularly in the *S 

case are lower than the corresponding 0 'Malley-Geltman states. 

The discrepancy in e 2 is 0.011 e.v. which, by comparison to (2.25) 

is seen to be a difference of 30$. The qualitative lowering of the 

levels as we have obtained is almost certainly correct. This is again 

because of the tremendously long range of the r 2 potential. A 

quantitative estimate of this potential can be obtained by observing 

when the second term in the asymptotic expansion (2.12) becomes 

comparable to the first. For any eigenfunction n, we define this 

value of r by R *. 

n 


R 


n 


a‘ 


2,1 


'n 


(2.24) 


In the 4 S case Ri — 100 ao where we have explicitly appended the 

unit of length ao (the Bohr radius). Having used a not ujjusual 

partial 

relative .wave type variational wave function with up to 25 terms, it 

\ 

is not surprising that O'Malley and Geltman could get an accurate 
lowest state whose range (25 a o ) is not very much more than an ordinary 
type bound state wave function. However, the range of the second 
state (100 ao) by a conventional variational wave function below the 
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n = 2 threshold at all. There can he little question that the value 
itself must he too high. 
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TABLE I 


Autoionization levels of H (in eV above the ground of hydrogen) 


State 

0 'Malley-Geltman 3 

This paper 


9-559 

9.559 


10.178 

10.1662 

10.2016 

3 S 

10.149 

10.149 


10.202 

10.2006 



10.2036 

*P 

10.178 

IO.I78 


10.203 

10.2029 


— 

10.2037 

3 P 

9.727 

9.727 


10.198 

IO.I875 


* 

10.2032 


a Our numbers are based upon the rydberg being taken as exactly 
13.605 eV. 


b From reference 1. 

c This column contains only the first three of the resonances 

for each state of the infinity of resonances predicted by Eq. ( 2 . 19 ). 



